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Collision Dynamics of Two Barchan Dunes Simulated Using a Simple Model 
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The collision processes of two crescentic dunes called barchans are systematically studied 
using a simple computer simulation model. The simulated processes, coalescence, ejection and 
reorganization, qualitatively correspond to those observed in a water tank experiment. Moreover 
we found the realized types of collision depend both on the mass ratio and on the lateral distance 
between barchans under initial conditions. A simple set of differential equations to describe the 
collision of one-dimensional (ID) dunes is introduced. 
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Many types of sand dunes are found in deserts, on 
the sea bottom and even on Mars. Dunes are formed 
by interactions between the flow of wind or water and 
sand. 1 The flow makes the shape of a dune by trans- 
porting sand particles. The wind transports sand parti- 
cles and forms the shape of dunes. The dune topogra- 
phy, in turn, acts as a boundary condition on the air 
flow. One of the most intensively studied types of sand 
dune is the crescentic dune called the barchan. A steadily 
blowing unidirectional wind generates barchans when the 
amount of available sand is insufficient for covering the 
entire bedrock. 2 Barchans usually migrate as a group, 
interacting with one another through collisions and in- 
terdune sand flow. 1, 3-11 Previous studies, however, have 
focused mainly on single barchans. 12-18 One of the rea- 
sons is that, because of the slow time scale of the sys- 
tem, it is hard to observe the whole process of the in- 
teraction dynamics of barchans in a desert. Recently, 
some underwater experiments have succeeded in creating 
barchans, 19-22 which has enabled real-time observation. 
In this letter, taking only dominant factors into account, 
we propose an effective model to simulate the coaxial and 
offset collisions of two barchans, which are observed in 
the water tank experiments. 22 

In the model, the dune field is divided into square 
cells. 3-5 Each cell is considered to represent an area of the 
ground which is sufficiently larger than the sand grains. 
A field variable h(x,y,t) which expresses the local sur- 
face height is assigned to each cell; t denotes the discrete 
time step and the spatial coordinates x and y denote the 
positions of the center of a cell in the flow and the lat- 
eral directions, respectively. The edge length of the cell 
is taken as a unit of length. In short, x, y and t are dis- 
crete variables while h{x,y,t) takes a continuous value. 
This model belongs to a class of simulation models called 
cell dynamics. 23 The motion of sand grains is realized 
by two processes: saltation and avalanche. Saltation is 
the transportation process of sand grains by flow. The 
saltation length and saltation mass are denoted L and 
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Fig. 1. (a) Relationships between height, length and width of 
barchans. The solid lines are the best linear fit. (b) Relation- 
ships between barchan velocity(V) and barchan height(H). The 
solid line is the best fit using the relation V = Q/(H + H c ) + V c , 
where H c and V c are phenomenological parameters. 24 



q, respectively. Saltation occurs only for cells on the up- 
wind face of dunes. In each time step of the simulation, 
the amount of sand q is shifted from a cell (x, y) to the 
leeward cell (x + L,y), which is the numerical expres- 
sion of a saltation process. Hence, the changes in height, 
h(x, y) — ► h(x, y) — q and h(x + L,y) — ► h(x + L,y) + q, 
take place at the taking-off cell and the landing cell, re- 
spectively. The saltation length L and the amount of 
transported sand q are modeled by the following rules, 



L = a + bh(x,y,t) — ch 2 (x,y,t), 
q = d. 



(1) 

(2) 

where a=1.0, b=I.O, c=0.01 and d=0.1 are phenomeno- 
logical parameters. In eq. (1), the second term shows that 
sand is transported farther away as the height of the sand 
surface is higher. The last term is introduced for L not to 
become too large. Note that (1) is used only in the range 
where L increases as a function of h(x,y,t). The salta- 
tion mass q is fixed at 0.1 for simplicity. In the avalanche 
process, on the other hand, the sand grains slide down 
along the locally steepest slope until the slope relaxes to 
be (or be lower than) the angle of repose which is set to 
be 34°. 

We examine if some basic features of barchan are real- 
ized by this simple model. 12, 14, 15 Starting from an initial 
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Fig. 2. Typical collisions of two barchans: (a) coalescence process simulated with offset parameter a = 0.0 and mass ratio j3 = 0.2, (b) 
ejection process simulated with offset parameter a = 0.0 and mass ratio fi = 0.7, (c) reorganization process simulated a = 0.3 and 
(i = 0.2, (d) offset collision performed in water tank which is 10 m long, 20 cm wide and 50 cm deep with water depth maintained at 
approximately 13 cm. 22 In (a), (b) and (c) the slopes at the angle of repose are painted in black. 



Gaussian sand pile, we measured the morphologic rela- 
tions after a steady barchan shape was reached. The lin- 
ear relationships between height, length (Fig. 1(a)) and 
width(Fig. 1(b)) are shown. Next, the roughly inverse 
relationship between migration velocity and height (Fig. 
1(c)) was confirmed to hold. In these tests, quasi-periodic 
boundary conditions were used in which the total mass of 
sand flowing away from the downwind and lateral bound- 
aries was homogeneously re-injected from the upwind 
boundary. 

In order to simulate collision processes, two initial 
barchans are situated in a field. The longitudinal (i.e, 
wind directional) distance between their crests is set to 
d x =2L i} where L ; is the main body length of the ini- 
tial leeward barchan, while the lateral distance between 
them is set to d y =aWi, where Wi is the width of the same 
leeward barchan. The constant a is varied as a control 
parameter which we call the offset parameter. We fix the 
initial mass of the leeward barchan (M;=10) and define 
its ratio (j3 = M w /Mi) to the initial mass of the wind- 
ward barchan (M w ) as another control parameter. 

In the first simulation, we examine the coaxial col- 
lisions. For the sake of quantitative comparison with 
experiments, open boundary conditions are used. This 
means that there is no influx sand from the upwind 
boundary and the sand out of the numerical field from 
the downwind boundary is neglected. Figure 2(a) shows 
the simulated process for a = 0.0 and [3 = 0.2. After two 
barchans collide with each other, the windward barchan 



Table 1. Calculated results classified into coalescence (x) , ejec- 
tion (A) and reorganization (o) for various mass ratios and offset 



parameters. 

Offset Mass ratio (/?) 

parameter (a) 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

xxxxxAA 

0.1 x x o o o o o 

0.2 x o o o o o o 

0.3 o o o o o o o 

0.4 o o o o o o o 

0.5 o o o o o o o 



is absorbed into the leeward barchan, during which the 
upwind slope of the leeward barchan is temporally eroded 
and forms a dent. As time proceeds, the dent is filled up 
and the unified dune recovers the crescent shape of the 
barchan. The behavior of barchans which is similar to 
this simulation result was observed in water tank experi- 
ments. 22 Hereafter, we call such a process the coalescence 
of two barchans. Figure 2(b) shows the time evolution for 
a = 0.0 and (3 — 0.7. Similar to the case of coalescence, in 
the initial stage, a part of the windward barchan merges 
into the leeward barchan. However, the valley between 
two crests is kept imperfectly filled, thereafter, a small 
barchan is bled from the downwind slope. Although the 
small barchan firstly drags ridges from its horns con- 
necting to the windward barchan, the ridges are soon 
cut and finally steady shapes of windward and leeward 
barchans are attained. Hereafter, we call such a process 
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Fig. 4. Schematic diagram for collision process of ID barchan col- 
lision. The solid and dashed lines represent the profiles at time 
t and time t + At, respectively. 
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Fig. 3. Time evolution of heights H= (Hi , H2 , H v ) in longitudi- 
nal profile along solid line as indicated in inset, (a) For coaxial 
collision (Fig. 2(a)) the coalescence of two crests is seen at ap- 
proximately 9000 time steps, (b) For offset collision (Fig. 2(c)) 
the exchange of the heights of two crests is seen at approxi- 
mately 9200 time steps, which indicates the reorganization of 
two barchans. 



the ejection of a small barchan. Figure 2(c) shows the 
time evolution for a = 0.3 and (3 — 0.2. As the initial 
windward barchan approaches the leeward barchan, the 
wing of the leeward barchan is partially eroded and is 
pushed downward. After that it is separated from the 
deforming main body. The separated body drags a ridge 
extending from the main body for a while. Subsequently, 
it is completely isolated and recovers the crescent shape. 
Hereafter we call such a process the reorganization of two 
barchans. This process is qualitatively similar to what 
happened in the experimental tank (Fig. 2(d)). Note that 
the types of realized collision depend on the mass ratios 
of initial barchans and on their relative positions. 

In addition, systematic calculations arc performed 
with various combinations of a and /3. The obtained colli- 
sions are classified into three types, coalescence, ejection 
and reorganization (Table I). Note that the coalescence 
appears in the region of low values of offset parameter 
and small mass ratios, and the reorganization occurs at 
high values of offset parameter or large mass ratios. On 
the other hand, ejection occurs at low values of offset pa- 
rameter and large mass ratios, which corresponds to the 
"solitary-wave behavior" obtained by the previous sim- 
ulation. 9 This diagram should be compared with those 
of further subaqueous experiments and also to those of 
field observations. 

To explore the intrinsic mechanism underlying the dif- 
ference between the three types of collision, we focus on 
the dynamics of longitudinal profiles of barchans partic- 
ularly on the heights of the windward crest (Hi ) , the lee- 
ward crest (H2) and the bottom of the valley (H v ). Fig- 



ure 3(a) shows the time evolution of H = (i/i, i/2, H v ) in 
the case of coalescence. H v reaches H\ at approximately 
9000 time steps. Figure 3(b) shows the time evolution 
of H in the case of reorganization. The exchange of the 
heights of two barchans is seen at approximately 9200 
time steps, whereas the bottom height of the valley is 
kept below them throughout the process. The time evo- 
lution in the case of ejection is qualitatively the same as 
that in the case of reorganization. In short, ejection is a 
type of reorganization. We discuss only coalescence and 
reorganization hereafter. 

The dynamics in each profile of Fig. 3 is not self- 
complete because of the lateral sand flow. Nevertheless, 
the dynamics realized in each profile is roughly inferred 
from the initial condition, that is, (i) if the height ratio of 
the windward crest to the leeward crest is comparatively 
small under the initial condition, coalescence occurs, (ii) 
if the height ratio is comparatively large, the reorganiza- 
tion of two barchans occurs. In process (i), the smaller 
barchan on the windward side climbs the upward face of 
the larger barchan on the leeward side and is absorbed on 
the latter before its crest reaches the same height as the 
latter. On the other hand, in process (ii), the windward 
barchan climbs the larger leeward barchan and its crest 
becomes higher than that of the leeward barchan with- 
out absorption into the latter. Subsequently, the leeward 
barchan runs away from the windward barchan because 
of the inverse relationship between velocity and height. 
From the above discussion, the important factor deter- 
mining the type of collision is presumed to be the com- 
petition between the absorption and height exchange of 
two barchans. 

To attain a clear picture of the complex processes of 
collisions, we further simplify the model using a set of 
naive assumptions. The first assumption is the geomet- 
ric similarity between ID dunes of different sizes which 
means that constant angles of upwind and downwind 
slopes are maintained irrespective of their size. The sec- 
ond assumption is that the sand flux at the crests of 
dunes is described by the product of eqs. (1) and (2). 
The third assumption is the mass conservation of sand 
crossing over a crest, and that all the sand blown across 
a crest becomes trapped in the downwind face of the 
same dune without directly escaping to the leeward. The 
combination of the above assumptions leads to a simple 
conservation relationship between the eroded sand on the 
upwind faces and the accumulated sand on the downwind 
faces of each dune(Fig. 4). The sand mass conservation 
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Fig. 5. Time evolution of heights of windward dunc(dashed), lee- 
ward dunc(solid) and valley(dotted) obtained by numerical cal- 
culation. The dash-dotted line represents time when 5h reaches 
Sh m i„ = 0.5. (a) Under the initial conditions of /ii(0) = 1.02, 
^2(0) = 10.0 and h v {0) = 0.01, coalescence collision occurs, (b) 
Under the initial conditions of /ii(0) = 6.00, /i2(0) = 10.0 and 
h v (0) = 0.01, the reorganization occurs. 

in the windward dune is expressed as, 

(fei + ftj) A (h. + h^-jK + K) ^ _ 

P 2 AXul = p 2 = 

(3) 

where hi and h v are the heights of the windward dune 
and the bottom of the valley at time t, respectively, 
whereas h[ and h' v indicate the heights of the windward 
dune and the bottom of the valley at time t + At. Ax u i 
and Axdi denote the horizontal displacements of the up- 
wind surface and the downwind surface of the windward 
dune during At, qi = q{h\) denotes the sand mass cross- 
ing over the crest per unit time, and p denotes the area 
density of the bulk of sand. Also, the change in height of 
the windward crest for At is expressed as, 

h[ -hi = AAx dl -AAx ul = Ah u (4) 

where A is the geometrical constant of similar triangles 
constituting ID dunes, namely, the ratio of their heights 
to the lengths of their bases. Using eqs. (3) and (4) and 
taking the limits At — > 0, Ax u i — > 0, Axdi — > and 
Ah\ — > 0, the evolution equation of the crest height, 

dhi q\A 1 1 

at p hi — h v hi 

is obtained. Similarly the evolution equations of the 
heights of the leeward dune and the valley are written 

as, 

dh 2 g2A 1 _L 

dt ~ p ( h 2 h 2 -h v h [ ' 

dhy = giA 1 g 2 A 1 

dt p [ hi - hj p h 2 -h v '' K ' 

In addition to the above closed set of differential equa- 
tions, we set the minimum depth as 5h m in- Here, the 
valley depth Sh is defined as Sh — hi — h v and two dunes 
are considered to coalesce once Sh reaches Sh m i n - Figure 
5(a) shows a typical collision obtained by the calculation 
of eqs. (5), (6) and (7) with 5h m i n = 0.5. In this case, 
the coalescence of two dunes occurs at the time repre- 
sented by the dash-dotted line. h v reaches hi before the 
exchange of heights occurs and coalescence occurs. This 



process is similar to that indicated in Fig. 3(a). Figure 
5(b) shows another typical collision under the same ini- 
tial conditions as indicated in Fig. 5(a) except for the 
higher windward initial dune. The exchange of heights 
occurs before h v reaches hi and reorganization occurs. 
This process is similar to that indicated in Fig. 3(b). 
This indicates that the set of simple differential equa- 
tions, eqs. (5), (6) and (7), contains the intrinsic features 
of collision dynamics. 

To summarize, the coaxial and offset collisions of two 
barchans are simulated using a simple model that in- 
cludes only saltation and avalanche processes without 
taking complex wind flow into account. Also we intro- 
duced a set of differential equations, which effectively 
describe the collision processes of ID dunes. 
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